NASA-CR— 1 97867 


Date: 

Submitted to : 
Title: 


TECHNICAL REPORT 


/// 

£) C/T~ 



January 21, 1995 

National Aeronautics and Space Administration 



5 ^ 


ACTIVE CONTROL OF PANEL VIBRATIONS 
INDUCED BY A BOUNDARY LAYER FLOW 


Grant Number: NASA NAG-1-1175 

Period Covered: August 31, 1990 — October 31, 1991 

Submitted by: Pao-Liu Chow (Principal Investigator) 

Professor and Chairman 
Department of Mathematics 
Wayne State University 
Detroit, Michigan 48202 


(NASA— CR-197867) ACTIVE CONTROL OF 
PANEL VIBRATIONS INDUCED BY A 
BOUNDARY LAYER FLOW Technical 
Report, 31 Auq. 1990 - 31 Oct. 1994 
(Wayne State Univ.) 30 p 


N95-23 1 82 


Unci as 


63/02 0042742 



1. Introduction 

In recent years, active and passive control of sound and vibration in aeroelastic 
structures have received a great deal of attention due to many potential applications to 
aerospace and other industries. There exists a great deal of research work done in this 
area. Most recent advances in the control of sound and vibration can be found in the 
conference proceedings [1, 2]. In this report we will summarize our research findings 
supported by the NASA grant NAG-1-1175. 

The problems of active and passive control of sound and vibration has been in- 
vestigated by many researchers for a number of years. However, few of the articles 
are concerned with the sound and vibration with flow-structure interaction. Exper- 
imental and numerical studies on the coupling between panel vibration and acoustic 
radiation due to flow excitation have been done by Macstrcllo and his associates at 
NASA/Langley Research Center (see e.g. [3,4]). Since the coupled system of nonlinear 
partial differential equations is formidable, an analytical solution to the full problem 
seems impossible. For this reason, we have to simplify the problem to that of the nonlin- 
ear panel vibration induced by a uniform flow or a bondarv-layer flow with a given wall 
pressure distribution. Based on this simplified model, we have been able to consider the 
control and stabilization of the nonlinear panel vibration, which have not been treated 
satisfactorily by other authors. Although the sound radiation has not been included, 
the vibration suppression will clearly reduce the sound radiation power from the panel. 
The major research findings will be presented in the next three sections. In Section two 
we shall describe our results on the boundary control of nonlinear panel vibration, with 
or without flow excitation. Sections three and four are concerned with some analytical 
and numerical results in the optimal control of the linear and nonlinear panel vibrations, 
respectively, excited by the flow pressure fluctuations. Finally, in Section five, we draw 
some conclusions from our research findings. 

2. Boundary Control of Nonlinear Panel Vibration 

Consider a rectangular panel whose mid-plane is bounded by 0 < x < a and 0 < y < 
6. A spatially uniform air-flow passes over the panel with a time-dependent mean flow 
velocity £/(£), (sec Fig. 1). For a long span b >> 1, when the transverse deflection w is 
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uniform in the y-direction, a one dimensional structural model for the panel vibration 
is used, (see Fig. 2). At a high flow speed, the linear piston theory for the aerodynamic 
forces is assumed to be valid. Then the panel vibration is described by the following 
nonlinear integro-diflerential equation: 


a?+ 4 ^ + /M=o, 


( 2 . 1 ) 


dt 1 • I- >'/ - g x < 

where m is the mass density of the panel, P is the compressive in-plane load, D = 
Eh 3 / 12(1 — z/ 2 ), E,h and v denote the Young’s modulus, the panel thickness and the 
Poisson’s ratio, respectively. The additional tension due to the panel stretching is given 
by the integral: 

r , dw , 


fd fj'll } 

A = {Eh /2a) / {—) 2 dx, 


and the aerodynamic force f(w) can be expressed as 

n > ~ pu ‘ f dw i 

^ W} M, ' dx + U dt ■ ’ 


( 2 . 2 ) 


[2.3) 


where p is the fluid density and M c is the flow Mach number. Let I\ and V denote, 
respectively, the kinetic energy and the potential energy defined by 

r , dw , 


ZC(u?) = J {^-) 2 dx, 

i Jo Ot 


VM = \j/{D( 


d 2 


w 


1 


dw . 


dx 2 


> + < 5 * - 


so that the total energy is given by 


E t = K t + V t . 


(2.0 


(2.5) 


( 2 . 6 ) 


To release the compressive force P, we first apply a tensile force Q at x = a. Then the 
net compressive force becomes 


R = (P — Q — N) 

and the modified equation (2.1) takes the form 

d 2 w 


(2.7) 


m 


dt 2 


R?E + D d ' w 


dx 2 


dx* 


= /(«>)• 


(2.8) 
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If we multiply Eq. (2.7) by w = and integrate the resulting equation over 0 < x < a, 
it yields the following energy equation: 


*<»> -£*<«> - <!><£)(£> -<*>£-*£)*} B 


dx 2 /v dx ' v dx 3 dx 
° dw , 


- £fWwdx-p£(-£)*dx. 
Since the left end is clamped, we have 


(2.9) 


3m (Pm 

— — 0 and — r = 0 at x = 0. 

OX OX 6 


( 2 . 10 ) 


At the right end, x = a, in addition to Q, we apply a twist moment M and a vertical 
point-force F so that the boundary condition takes the form (see Fig. 3) 




( 2 . 11 ) 


The objective is to choose the control {Q, M , F ) properly to stabilize the system. 

By introducing the scalings x = x/a. i — t(D j ma 4 ) 1 ' 2 and w — wjh etc., the 
initial-boundary value problem for the boundary control can be written in the following 
dimensionless form: 

d 2 w ~d 2 w d*w 


di 2 


- R- 


= f( w), 0 < x < 1, 


dx 2 dx 4 

dw . . n 

u; = u'o and — r = uq at t = 0, 

oS‘ 

w — 0 and -gT = 0 at x = 0, 

d 2 w i> j b d ™\ b * * 1 

— = -11/ and (—-R-) = F at x = 1, 


( 2 . 12 ) 


where 


with P = P/a 2 ,Q = Q/a 2 and 


R = (N + Q — P) 


r 1 dw 


N = Q / (??) 2 dx, a = Eh 3 /2a 2 , 
Jo dx 


* . -dw dw 

/M = -(% + i'ar)- 


The remaining non-dimensional quantities are given by: /? = pa 3 U 2 / DM C , 

7 = pa 2 U/Mc{DM) 1 /2 , u\ = uq//i. i = 0, 1; A/ = AIa 2 /Dh and F = Fa 3 /Dh. 
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Next we consider the case of free vibration without flow excitation, where the com- 


pressive force P is static. For simplicity, from now on, we shall drop the caret symbol 
over the dimensionless variables. For no flow, we have / = 0 so that Eq. (2.12) yields 


d 2 w d 2 tw 

IF - R 1F t a* = °' 0<x<1 ’ 

ow 

w = Wo and = tui at t = U, 

dS,‘ 

w = 0 and = 0 at x = 0, 

d 2 W nr J /<P' W r>d W \ r? *. 

— =-A/ and (--R—) = F at x 


= 1, 


where 


(2.13) 


R = AP + N with A P = (Q- P), 



The corresponding kinetic energy A and potential energy V are 


I< 

V 


1 ri 

I\ t {w) = - u 
L Jo 

aH = l/ o « 


W 2 dXn 


It can be shown that, in order to have a positive-definite V, we require A P + 1 > 0 or 
Q must be chosen so that 


[P — Q) < (1 — 6) < 1 for any S > 0. (2-14) 

For P < 1, we simply taken Q = 0, or no tensile control is needed in this case. One 
notes that Per = 1 is the critical buckling load for a simply supported panel. 

In view of the boundary conditions in Eq. (2.13) and the fact / = P = 0, the energy 
equation (9) yields 

E t (w) = Fvi{t) — (2.15) 

where t)j(f) = w(t, 1) and 0,(f) = f &(t, 1). For the energy E to decay, we will choose 
a feedback control pair of the form F = g(v i) and M — h(0 1 ) such that E t < 0. An 
obvious choice is to assume g and h being linear, 


F = -fivi 


and 


M = vOx, 


(2.16) 
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where fi and v are some positive constants. They can be regarded as the damping 
coefficients against the right end’s transverse and rotational motions, correspondingly. 
A substitution of (2.16) into (2.15) gives 

E t = -pvl-v0l<O (2.17) 

This shows that, if the sub-critical load condition (2.14) is satisfied, the linear control 
(2.15) will results in a decay in energy. In fact the rate of decay is exponential. To 
verify this fact, we need to introduce a perturbed energy Ef as follows 

E e t (xv) = E t (w) + eG t (w), (2.18) 


where 0 < £ < 1 and 


= i‘ I( 


First, it is easy to show that 


So that Eq. (2. IS) leads to 


|G,(u’)| < E t (w) 


(1 -s)E t < E‘ t < (1 +e)E u 


(2.19) 


( 2 . 20 ) 


or, E t and E f s , as far as the exponential decay is concerned, are equivalent. Next, by 
differentiating Eq. (2.19), invoking Eq. (2.13) and some mathematical inequalities, we 
can derive the following inequality: 


E‘ = E t + £G t 


1 


< —eEt - [fi - ^(1 + n 2 )e]v\ - v{\ - ve)6\. 

If we choose e in the range 0 < £ < £o < 1 with £q = min{l,2/x/(l + /x 2 ),l/i/}, and 
make use of (2.20), it follows from the above inequality that 


El<~ 


which implies 


(1 + e) 

Et < Et,e~ xt 


Et, 
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or, by noting (2.20) 


E t < ( ■ ^ £ )^o e A< i A = e/(l + e). 

Thus we have shown that the free vibration can be exponentially stabilized by a small 
boundary damping together with a tensile control Q , which is necessary only if the 
compressive force P exceeds a critical level specified by (2.14). 


Now we consider the case of flow-induced vibration. Dropping the caret symbol over 
the dimensionless quantities again, under flow excitation the controlled problem (2.12) 
reads as follows: 


d 2 w d 2 w d 4 w dw dw 
, dw 

w — w o and — - = uq at t = U. 
n ot 

OIL' 

w = 0 and — — = 0 at x = 0, 

= -M L " 


( 2 . 21 ) 


dx 2 


w _ Ow . „ 

a?-- R T I ) = F at I = 1> 


where R = (N + Q — P). ,3 and 7 are time-dependent. By making use of Eq. (2.21) and 
integrating by parts, we can obtain the energy equation 

E, = AP f'&fdx -m f'wAdx- 7 (() f'&’dx (2.22) 

J 0 dx Jo ox Jo 

+ Fi'iit)- 


where A P = (Q — P). 

First we choose the control law (2.16) and (2.14) for (P. M) as in the previous 
case. Then, by applying some integral inequalities, Eq. (2.22) yields the following 
inequality: 

E, < (A> + /J»/2 7 ) J f(^) J *r-| J fd , dx 
- nvl - i/O*, 


which is negative if 


Thus it is sufficient to set 


A P + (3 2 j 2 7 < 0. 


Q < (P — /3 2 /27 — 6 ) for any 6 > 0. 


(2.23) 
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Then 


Et < " f w 2 dx - 6 f (nr~) 2 dx 
2 Jo J o ox 

- (IV 2 - vO\. (2.24) 

Let E\ = E t + eG t as before. In this case, one can verify that 

P, < [wHx-S J\^fdx 

~ it* ~ ^(1 + ^ 2 )eR 2 ~ ^(1 - (2-25) 

< -sE t , 

provided that 0 < £ < £i < 1 and 

£ X = min{1, 2/i/(l + // 2 ), l/^,max[2/?(0/7(0]}- ( 2 -26) 

As in the free-vibration case, by noting (2.20), we conclude that 

E t < ( | — )£’oe~' u , with A = e/(l+e), 

which shows the exponential decay of energy by the given boundary control. 

We remark that 

1. By a more careful examination of the conditions (2.14), (2.16) and (2.26) for stability, 
they imply that the flow velocity can only have weak fluctuation which tends to zero 
quickly as t — ► oo. 

2. In lieu of (2.27), the exponential stability holds if 

max{|A/ 3 | + /? 2 /2 7 } = £ 2 < 1. (2.27) 

This condition means that both the rate of change of the difference (Q — P) and the 
flow velocity are small. Since the governing equation holds for a high speed flow, the 
condition is unrealistic. 

3. In general, the boundary control as discussed here may not be sufficient to stabilize 
the system. A more active control is required for this purpose. 
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To illustrate the analytical results, we consider the controlled system (2.21) in the 
special form: 


where 


d 2 w dw dw 

~dF ~ R ~d^ + ^ Ih + 1 ~dt ~ 

iy(0,x) = “(1 - cos 2ttx) and — (0,x) = 0, 
dw 

w — 0 and — — = 0 at x = 0, 

Oil) 

d 2 w d~w 1 (Pw dw 

= and d^- R d! = "aT at I = 1 ’ 




(2.28) 


(2.29) 


The associated energy of the system is given by 

1 /•! #u\ 0 1„.3 u! >0 ,d 2 i 


, 1 r l r/ oiL \ 0 1 .o w.r,. , 

b <°- 5JW + 2* ( & ) ' + ( &r )) *- 


(2.30) 


To show the effect of boundary damping on the panel vibration, the system energy 
was computed based on a modal expansion. Since the boundary condition in (2.23) 
is non-standard, the eigen-function expansion for the static problem cannot be used. 
Instead we approximate the deflection w(t,x) by a truncated Fourier series as follows: 


1 


A r 


w(t,x) ~ -a 0 {t) + ^{a n (0 cos 2n~x + b n {t ) sin 2 n~x}, 

~ n=l 


(2.31) 


where N > 1 is a fixed integer. When the above series is substituted into Eqs. (2.28), 
we obtain a coupled system of (2 N +1) ordinary differential equations of the form: 


a n (0 = /jv(<*o, a u ■ ■ ■ , a A'i bi,..., fc-v), «n(0) = Q„,n = 0, 1, . . . , A, 

^n(0 = 9N{ a 0i a l ? • • • i a N‘i b\ , • • ■ , ^n(O) = Pm Tl = 1,2..., A, 

The above system was solved numerically by the 4th-order Runge-Kutta method for 
N = 2, 3, 4, 5. We found that for N — 5, the modal amplitudes a n and b n are numerically 
negligible for n > 3. Thus all the results to be shown were obtained by truncating the 
series (2.25) at N = 3. The corresponding system energy E(t) given by (2.24) was 
evaluated at i/ 0 = 0-33 and various other parameter values. The numerical results were 
summarized and displayed in Figs. 4-10, which exhibit the evolutions of the vibrational 
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energy levels with and without boundary damping (control). In the case of free vibration 
(/? = 7 = 0), Fig. 4 shows that, when /x == v = 0 (without control) the energy is at 
the constant level of 130 and, when /x = u = 1.5 (with control), it does indeed decay 
exponentially, in agreement with the theoretical prediction. For the ease of visualization, 
in the sequent figures, the results, with flow excitations, for the uncontrolled and the 
controlled cases were plotted separately. In Fig. 5 with flow parameters /3 = 50 and 
7 = 0 (without aerodynamic damping), the energy level oscillates periodically when 
there is no boundary damping. By contrast, given the control parameters p = 2.5 and 
v = 1.5. Fig. 6 shows the exponential decay of the system energy. Corresponding to 
Figs. 5 and 6 , similar results are displayed in Figs. 7 and 8 when the flow damping 
parameter value was changed slightly from 7 = 0 to 7 = 0.05. The most drastic effect 
of the boundary damping is shown by Figs. 9 and 10 where the flow speed is high 
(j3 = 500) but the flow damping 7 = 0 . Without control, the energy level osciallates 
very rapidly, as seen from Fig. 9. This highly oscillatory state of energy can be reduced 
to an over-damped state, as depicted in Fig. 10, by introducing a boundary damping 
with control parameters p = 14 and v = 6 . 

3. Optimal Control of Linear Panel Vibration 

Let us consider a viscous flow past over the elastic panel. The flow is governed bv 
the well know Navier-Stokes equation: 

^ + (u-V)u = --Vp + */Au, (3.1) 

at p 

where the notations are standard. For a slightly compressible flow, the continuity 
equation reads 

j- t + V . (pu) = 0. (3.2) 

The panel is regarded as an elastic plate with thickness h. Young’s modulus E and 
Poisson’s ratio 7 . Under a uniform tension with T > 0 (or compression with T < 0) 
and the fluctuating wall pressure, the vertical displacement £ of the plate satisfies the 
following equation: 



where p w is the plate density, p w the wall pressure fluctuation, and q is the applied force 
as the active control. The constant D is the stiffness of the plate defined by 


D = 


Eh 3 


(3.4) 


12(1 -7 2 ) 

According to the boundary-layer theory, given an upstream velocity field U , the 

flow near the plate can be determined by the Prandtl’s approximation. In particular, if 

dp 

the panel is located on the x — y plane, the pressure gradient — across the boundary- 
layer is nearly constant, where x = (x,p, 2 ). Suppose that the mean flow outside the 
boundary-layer is parallel to the plate so that U = (^oo,0,0) + U\{x,t ), where U\ is a 
small perturbation. To derive the equations for acoustic quantities and p \ , we let 


u = Uq + 1 * 1 , p = po + Pi and p = p 0 + pi 


(3.5) 


where u 0 ,p 0 anc l Po are flow variables associated with the mean flow. As in the stability 
analysis, we introduce a parallel flow approximation. Then, in view of (3.5), one obtains 
the acoustic equations from (3.1) and (3.2) by linearization: 

^7 + (fio-V)u 1 = --V Pl , (3.6) 

at Po 


■J— + V • (p Q ui + P\Uq) — 0. (3-0 

For an isentropic flow, p\ and p\ are related by 

Pi =p>/c ! , (3.8) 

where c is the speed of sound for the unperturbed flow. Aside from a static displacement, 
the vibration of the panel is described by the perturbation w of equation (3.3) as follows: 

pj : = TAw - DA 2 w + f + q(x,t), (3.9) 

where f = p w is the fluctuating part of the wall pressure excited by the unsteady 
boundary-layer flow. The coupling of the acoustic equations (3.6) and (3.7), and the 
plate equation is through the boundary conditions. For the plate equation (3.9), since 
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it is simply supported by a periodic structure, we need only to analyze the problem 
over a fundamental domain 0<x<a,0<j/<6 and impose the boundary conditons: 


w(x,y,t) = 0 at x — 0,a; y = 0, b. 


(3.10) 


dp 

Since the pressure gradient — = 0 across the boundary, the wall pressure p w can be 
determined from the perturbed potential flow field Ui through an approximate Euler’s 
equation, that is, 


Pw = F{U X ). 


(3.11) 


To counter this excitation, a control force q(x,y,t) was introduced in (3.9). The objec- 
tive of the active control is to minimize the average vibrational energy and the control 
cost: 

J{q) = 2 tJo Jd {q{ JF )2 + 31 ^ + 7 ' Vu ’ i2 + k< ?} dtdxdy ' ( 3 * 12 ) 

where the time T may be infinite, D is the basic domain {0<x<a,0<?/< 6}; ct,3, 7 
and k are positive constants. In the language of the optimal control of a distributed 
parameter system, the equation (3.9) is known as the equation of state and J(q ) defined 
by (3.12), the objective or cost functional. Here the physical problem of vibrational 
control reduces to an optimization problem: Given the wall pressure excitation p w , 
find an optimal control q*(x,y,t) from a certain admissible class Q of functions which 
minimizes the objective functional J(q ), that is, 


J(q*) = min{ J(q), q in Q}. 


(3.13) 


To obtain an analytical solution, we consider the case of simply supported boundary 
conditions: 

w(x,y, t) — 0 at x = 0, a; y = 0, f>, 


d^w d^w 

-Q^{x,y,t) = 0 at x = 0,a and ~^(x,y,t) = 0 at y = 0, b. 


(3.14) 


The initial conditions are given by 


dw 


w{x,y, 0) = W 0 (x,y), ~(x,y, 0) = uq(x,y). 


(3.15) 
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It is well known that the set of functions 


/ \ • rrz 7i r , ^ / o i c? \ 

^mn(x,2/) = 2sin— xsm— j/, m,n = 1,2,- •* (J.lbj 

are orthogonal eigenfunctions associated with the plate equation (3.9) and the corre- 
sponding eigenvalues are 


(3.17) 


In terms of the above eigenfunctions, we can expand the displacement te, the wall 
pressure / and the control as follows: 

CO 

^ ^m t n (Or^mn i !/ ) ? (3.18) 

m,n= 1 


CO 

= ^ y /m,n (O^ron. 2/ 7 

m ,n= 1 


and 


oc 

^ ] 9m,n (O^mn I/)? 
m,n= 1 

where the coefficients te mn etc. are computed by 


(3.19) 


(3.20) 


Vu\ 


in — — / / ^'(^b lb 0 ( r ? mn(‘ Z b 

Jo Jo 


and so on. A substitution of the expansions (3.18)— (3.20) into the equations (3.9), 
(3.15) and (3.12) yields the following uncoupled system of equations: 


Pw^mn “1" ^mn^mn fm n ( 0 9mn(0 5 

(0 ) = ^0,mn) ^mn(O) ^l,mn 


(3.21) 


and 


00 

J(q) = E , (3.22) 

m,n = l 

where 


1 T 

T j Q { a ^mn( 0 + / i ">n^L(0 + ^'ZmnCO}^ 


(3.23) 



for m,n = 1,2,*-*. Since the modes are uncoupled, if the cost J mn for each mode is 
minimized, so does the total cost J. 

For a given (m,n)-mode, dropping all the subscripts, we are led to consider the so- 
called “linear regulator” problem in optimal control: Find the control q in the equation 

j pw + A w = f(t) + q(t) , (3.24) 

\ w( 0 ) = w Q , w(Q) = W X , 

which minimizes 

J{q) = f ( Q ^ 2 + / m ’ 2 + kq 2 }dt , (3.25) 

11 JO 

where p mn is given as in A mn with D and T replaced by (3 and 7, respectively. By the 
method of adjoint state, for the cost to be minimal, the state w and its adjoint v 
must satisfy the optimality system: 

f pw + A w = f(t) + jfau’ — v ), (3.26) 

\ io(0) = ic(0), w(0) = u’i, 

and 

f pv -f Av = (aA + p)w, /g 9 -\ 

\ v{T) = v{T) = 0. V 

The optimal control q* is given by q = A (ate — u). One notes that, due to the coupling 
between w and 1?, the above system (3.26)— (3.27) is a two-point boundary- value problem. 
Numerically it can be solved by the shooting method. Some numerical results for the 
original modal equations (3.21) and (3.23) have been obtained. 

For example, we choose a — 47 t, b = 7T, Wq = =0 and T = 4, and set 

C0S ^ m2 + n2 ) 1/2f ’ n = 1, 2 , . . . . 

All the parameters are taken to be one except for (3, which is zero. The maximal 
amplitude of vibration under a optimal control has been computed and some results, 
corresponding to 4 modes (m + n = 4), are shown in Fig. 11 to Fig. 13. In the above 
figures, the solid curves represent the controlled amplitudes, which are in contrast with 
the uncontrolled ones. It is seen that the control is very effective in reducing the 
vibration amplitudes. For an independent interest, the controlled mode shape at t = 4 
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is plotted as shown in Fig. 14. 

4. Optimal Control of Nonlinear Panel Vibration 

Similar to Eq. (2.1), under a tension T and without flow excitation, a simplified 
model of the control of nonlinear penal vibration, in lieu of Eq. (3.9), is given by 

(4.1) 


= {T + /V)ttt - “I + /(*»*) + 9( x > 0 .° < x < G - 


d t 2 


dx 2 


dx 4 


with w(x,0) = g(x), and |^(x,0) = h(x), where, as before, 


.. Eh f a Ow 2 
N = — (—Ydx. 
2 a Jo ' dx 


dw . 


For a simply supported beam, we have 

d 2 w 


w = 0 and 


dx 2 


= 0 at x = 0. a 


while, for a clumped beam, t lie following holds 

c)iv 

w = 0 and — — = 0 at x = 0. a. 
dx 


(4.2) 


(4.3) 


(4.4) 


Other boundary condition are possible. In the subsequent analysis, we consider the 
simply supported case (4.3) only. The objective is to choose an optimal control q*(z,t) 
which minimizes the following cost function: 


j(q) = r Aafe 2 + 0(j%) + “ 2 + <>*■ 

2 In Jo Jo of CJZ 


(4.5) 


where a,/? and k are positive constants. 

The necessary condition for minimum is the vanishing of the variation SJ of J, or 
CT , x 1 f T f a r ,di u„d8w\ , a/ d 2 w^ ( d 2 8w y , 

SJ{q) = tL L {a ^ )( w )+ p { w )i ~w )+wSw+kqSq] (4 - 6) 

dtclx — 0, 

where 6w and 6q denote the variations of w and q , respectively. Note that, by taking 
the variaiton of equations (4.1) and (4.3), we can relate 8w to 8q by the variational 
equation: 


... d 2 8w d 2 8w d 4 8w ,9 2 u; 

M 8w = p-xrr = (T + N t )-^r + D-^- - ( 8N t )-^7 = Sq, 


dt 2 


dx 2 


dx 4 


dx 2 


(4.7) 


B 


8w = — Sw = 0 at t = 0, 
dt 


)W 


d 2 

dx 2 


8w = 0 at x — 0, a, 
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where 


r a .dw..dSw. , 

SN ■ = -< L ( T x ){ -SI )ix ' 

Eh 


(4.8) 


with 7 = 

a 

Upon eliminating 6 q from Eq. (4.6) by Eq. (4.7) and noting Eq. (4.1), it yields 

SJ=l [ T f {(-<fe + 0^ + w)8w + kqM5w}dt dx = 0, 

T Jo Jo at 2 ox 4 

which implies, after integrating by parts several times, that the optimal control q must 
satisfy the following equation. 

A , nA , r 


d q _ _ v £/ <7 _cr 4 o criu Gritf 

A/ * 9 = _ + iV fe + + 7 J 0 q j^- 


dx 


(4.9) 


1. d 2 w > . 

= Tl a -gF ~ ~ W) = ' 


5.i 


and the terminal - boundary condition: 


q 

q 


0 and ^ + 4^ = 0 at t - T, 
at pk at 


: A 

dx 2 


= 0 at x = 0 , a. 


(4.10) 


The equations ( 4 . 1 ) and (4.9) together with the conditions (4.3) or (4.4) and (4.10) form 
the optimality system. A peculiar feature of the system is that the state w is coupled 
to the control q which satisfies a terminal condition. This has caused some difficulty in 
computing the optimal solution. 

To solve the system, we proceed by applying the Galerkin method of approximation. 
For the simply supported case (4.3), we may use the following set of admissible function 

yj n (x)\4 sin— ,n = 1,2,... (4.11) 

V a a 


as an orthonormal basis. Now we expand w and q as follows: 

oo 

W(x,t ) = U>n(t)<Pn{x), (4-12) 

n = l 

oc 

0 = ?n(0 { r ? n (^) • 

n= 1 
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Substituting the series (4.12) into equs. (4.1) and (4.9) yields 

mw n (t) + [(T + Q) + D\l]\lw n = f n (t) + q n (t), (4.13) 

m 9n(0 + [(T + Q) + D 2 \ 2 n \\^q n 

OO ^ 

+ 7( ~ ~ ( CXU'n ~~ 0\ n W n ~ U>' n ), Tl = 1,2,..., 

m = l 

^n(O) = g n ,w n (0) = h n , (4.14) 

r\ 

<ln(T) = 0,).(r) + - 1 i'„(TJ = 0, 

where = j t w n) J n = | 0 a f l ^ n dx i etc; A„ = ^ ere! 

OO 

Q = £ 

m = 1 

By a modal truncation, the above infinite coupled system can be reduced to a finite 
system as follows: 

rml\V,n + [{T + Qx) + D\l]\lwi\\n = fn + qn, (4.15) 

N 

m ^V,n 4* [(T + Qn) 4 tfiV.n + ^l^Xm?A>) 

m = l 

= “[azD.v.n = (1 + /9A^)u?.v ?n ]? n = 1, 2, . . . , iV, 

which are subject to the conditions as in (14), where Qn is iV-term truncation of Q. 

The truncated system (4.15) can be solved numerically. The numerical solution of 
the truncated problem poses two technical difficulties: the large scale in computation 
and resulting two-point boundary values given in time domain. As a result, we can only 
handle a small set of modal equations. To be computationally efficient, we adopted the 
so called shooting method. The method consists of solving the boundary value problem 
as an initial value problem by assigning the missing initial data and then adjusting the 
data by interations until the end point conditions are met. The interation procedure is 
based on a fixed point algorithm in locating the zeros of a function. 

The numerical computation has been carried out in a SUN workstation. The results 
confirmed the theoretical prediction that the feedback control can drastically reduce the 
panel fluttering and it is more effective in suppressing the lower frequence vibrations. As 
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an illustration, by suitably choosing the parameter values, some results for a three mode 
calculation are shown in Figs. 11-13 under a time-harmonic wall pressure fluctuation. 
For the given numerical results, we set N = 3 ,/,(£) = (sinf/2 + cos£/2)/3i, and the 
parameter values: m — 1, a = 0.1,/? = 1,7 = 0.01, D = 0.02, k = 1. The domain 
of computation is given by 0 < x < 4 and 0 < t < 4. The amplitudes of the first 
three modes Z U Z 2 and Z 3 are plotted in Figs. 11-13, where the solid curves depict 
the amplitudes for the controlled system in contrast with the chain- like curves for the 
uncontrolled case. The reduction in the modal vibration amplitudes is very remarkable. 
Consequently, if we consider the sound generated by the vibrating panel, the sound 
radiation power or the aero-acoustic noise level will be reduced drastically by the active 
control. 

5. Conclusions 

Several problems in the control of panel vibration due to the flow excitation were 
studied. The main research results are presented in the previous three sections. The 
problems were solved analytically and numerically. Based on these results, we can draw 
the following conclusions: 

In the boundary control of nonlinear panel vibration, by means of the energy method 
and some mathematical inequalities, the boundary stabilization of a vibrating nonlinear 
elastic panel was studied. The panel is clamped at one edge and free to vibrate at the 
other edge. In general, the panel is subject to a compresssive in-plane loading combined 
with an aerodynamic forcing. Without any control, the panel would flutter due to flow 
induced instability. To stabilize the panel, a boundary control was introduced as the 
combination of a bending moment, a vertical point force and a tensile force applied to the 
free edge. Two cases, corresponding to the absence and the presence of an aerodynamic 
loading, were treated separately. For no flow, this is the case of free vibration. Even 
though the energy of the uncontrolled system is conserved, with initial disturbance, the 
system may buckle or sustain a persistent large-amplitude oscillation. To render the 
energy an exponential decay, it was found sufficient to apply a tensile force, if necessary, 
to reduce the net force to a subcritical level and, at the same time, to introduce a 
boundary damping. The damping mechanism consists of a pair of frictional force and 
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torgue, which are linearly proportional to the transverse and the angular velocities 
of the right edge, respectively. Therefore, if the compressive force is subcritical, the 
passive control in the form of a boundary damping suffices to stabilize the system. In an 
analogous situation, the result seems to be in agreement with the experimental evidence 
that a boundary damping is effective in suppressing the panel vibration [6]. In the 
presence of unsteady flow and compressive force, the panel is subject to flow excitation. 
If the flow velocity is oscillatory and decays rapidly, it is possible to stabilize the panel by 
applying a time varying tensile force Q{t) together with a boundary damping as before. 
But the control force Q must follow the flow fluctuation closely. For this reason, such 
scheme is not robustic. However, for slowly varying compressive force and, at the same 
time, small flow parameter 5 2 /7, the system can be stabilized as in the free-vibration 
case. It is believed that, when the flow parameter is large, the boundary control may 
not be sufficient for stabilization and a stronger mode of control, such as a distributed 
control, will be required for this purpose. The preliminary results of this results were 
summarized in a paper which was published in Recent Advances in Active Control of 
Sound and Vibration, Volume 2, [7]. The full paper [S] containing the analytical and 
numerical results was accepted by the Journal of Sound for publication and Vibration 
and will appear shortly. 

In the control of both linear and nonlinear panel vibrations, we investigated some 
optimal control problems with and without flow excitations. For the linear problems, 
the vibration control of a simply - supported rectangular plate was treated. The con- 
trol objective is to minimize the objective function, which is the sum of the vibrational 
energy and the cost of control. For the nonlinear control problem, a simplified panel 
equation similar to the boundary control problem was used. By the optimal control 
theory for distributed parameter systems, the problems can be solved by deriving the 
optimality equations for the adjoint states. These equations coupled with the dynamical 
equations for the panels must be solved to yield the optimal control forces. For approx- 
imate solutions, we adopted the Galerkin method or the modal expansion to reduce 
the governing partial differential equations to a finite system of ordinary differential 
equations by truncation. For the linear elastic plate, the control of each modal ampli- 
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tude reduced to the so-called linear regular problem which can be solved with relative 
ease. For the nonlinear control problem, the modal equations and their adjoint-state 
equations are nonlinearly coupled and the solutions are difficult to obtain. However 
numerical solutions were carried out. The numerical results show clearly that the con- 
trol is very effective in suppressing the panel vibration and the sound generation by the 
elastic panel. The results will be summarized in one or two papers to be submitted for 
publication. 
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Figure 2: Schematic Diagram of Panel Vibration with Flow Excitation 
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Figure 4: flow parameters: 3 = 0 , 7 = 0, uncontrolled and controlled with 


li = 1.5, v = 1.5 
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Figure 11: Amplitudes of controlled and uncontrolled vibration at the center of the plate 
with m = 2, n = 2 and T = 4. 
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Figure 12. Amplitudes of controlled and uncontrolled vibration at the center of the plate 
with m = 3,n = l and T = 4. 
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Figure 13: Amplitudes of controlled and uncontrolled vibration at the center of plate with 
m = l,n = 3 and T = 4. 
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Figure 15: The controlled and uncontrolled amplitudes of the first mode at the center with 
m = 3 and T = 4. 
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Figure 16. The controlled and uncontrolled amplitudes of the second mode at the center 
with m = 3 and T = 4. 
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Figure 17: The controlled and uncontrolled amplitudes of the third mode at the center with 
m = 3 and T = 4. 




